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Abstract 

A Hamiltonian -ff(r) applicable to cuprate HTS, with a doping 
dependent pairing interaction T{x) ~ V{x) + U{x), is linked to a 
Cu3d-02p state probability inodel(SPM). A consequence of doping 
induced electron hopping, the SPM mandates that plaquettes with 
net charge and spin form in the CuO plane, establishing an effec- 
tive spin-singlet exchange interaction U{x). The U{x) is determined 
from a set of probability functions that characterize the occupation 
of the single particle states. An exact treatment of the average static 
fluctuation part of H shows that diagonal matrix elements f/kk < 
produce very effective pairing, with significant deviation from the 
mean field approximation, which also depends on a phonon-mediated 
interaction V. This deviation is primarily responsible for the diverse 
set of HTS properties. The SC phase transition boundary Tc{x), the 
SC gap A (a;), and the pseudogap Apg{x) are fundamentally related. 
Predictions are in excellent agreement with experiment, and a new 
class of HTS materials is proposed. Large static fiuctuation results 
in extreme HTS and quantum criticality. 
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74.72.-h 
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1 INTRODUCTION 

The 1986 discovery of non-elemental high Tc supercon- 
ductor(HTS) cuprates[l introduced a class of antiferro- 
magnetic Mott insulator ceramics that exhibit an exotic 
array of seemingly disparate superconducting(SC) and 
normal state(NS) properties. [H [3] Despite the plethora 
of experimental data indicating several microscopic in- 
teractions, there is incomplete agreement on which in- 
teractions are essential and how they combine to pro- 
duce HTS.[4j No current Hamiltonian leads to the di- 
verse HTS properties observed in cuprates. This paper 
presents a Hamiltonian, conjoined with a state probabil- 
ity model(SPM), that predicts the observed properties 
listed below and elucidates a key ingredient for pairing 
glue. It is shown that the pseudogap results from doping 
induced charge-spin fluctuation, and a new class of HTS 
materials with an anti-pseudogap is suggested. 

The disparity between elemental low Tc < lOK su- 
perconductor(LTS) and HTS properties is remarkable. 
The doped cuprates are extreme type-H superconduc- 
tors, exhibiting strong magnetic field induced quantum 
fiuctuations, with the interpretation that the state of 
these materials may lie in close proximity to a quan- 
tum critical state. [5]!^ Denoting x as the hole (electron) 
doping concentration, nine diverse cuprate HTS proper- 
ties are: order of magnitude increased a) average en- 
ergy gap amplitude A(x,T), and b) Tc{x), which is 
maximum at optimal doping Xop, can exceed lOOK, 
with Tc/Tp ^ 10^1 - 10^2. c) large shape ratio 
A{xop,0)/kBTc{xop) ^ 2, the maximum BCS LTS 
value, d) inverted parabolic (dome) shaped phase bound- 
ary Tc{x), e) magnetic field penetration with A^^(T) 
linearly decreasing with T for T << Tc, f) x dependent 
evolution of the isotope effect, g) a large condensation 
energy A^{x), and an anomalously large discontinuity 
in the electron specific heat at Tc{xop), h) an SC gap 
Ak with d-wave symmetry, contrasting the s-wave gap 
in LTS materials, i) a tentative first order phase transi- 
tion, with a concomitant quantum critical point. 

In the normal state an anomalous property is the 
suppression of the density of electronic states referred 
to as the pseudogap. The pseudogap Apg{x), estab- 
lished by a number of spectroscopic probes, summa- 
rized in Ref. [7], deceases linearly in the underdoped 
domain from a maximum value Apg{x 0.05) until it 



merges on the overdoped side x > Xop with the SC gap 
A{x). Some researchers refer to Apg as the SC gap, 
but there is considerable experimental evidence that the 
distinct lower energy SC gap exists, with a very dif- 
ferent a; < Xop dependence. [3 |H] Current theories ei- 
ther promote the pseudogap state as a pairing precur- 
sor of the SC state, or as a competitor due to unre- 
lated dynamical fluctuations. [U [10] Hence, the origin 
and the effect of the pseudogap on the SC state remains 
unresolved, pi [71 [TTl [5] The SPM gives new insight here. 

Cuprate HTS theories abound. Partially sup- 
ported by experiment, various theories and reviews 
thereof include electron-phonon interactions [12' with 
possible small polaron.|131 IT4] or Jahn- Teller polaron 
formation. |15) interband interactions. |16[[T7| exciton me- 
diated interactions, [HIIllllDl mi negative U-center 
pairing, [231 [H] bosonic electron- hole pairing. [23] spin- 
exciton,[25] spin-phonon interactions. [?f] fractals, [2H] 
quantum oscillations in fermi liquids. [29] quantum 
criticality,[3ni [SD 13 [321 [S] 2D strong electron-electron 
correlation with resonating valence bonds(RVB) and 
spin exchange, fluctuating spin exchange (overdoped 
regime), spinous, holons, various slave-particle tech- 
niques with extensive monte-carlo calculation, etc. [33l 
IHll35l[a[36l|37l[38l[3ll39l|40l|411[42] The RVB 
theory, 03] generally implemented using the Hubbard 
t — J Hamiltonian, [33] gives a reasonable perspective of 
the undoped antiferromagnetic charge transfer insulator 
phase, and as a potential HTS model, spin exchange has 
appeal since d-wave pairing is a natural consequence. [43] 
Although various HTS characteristics are explained by 
reasonable, but disjointed, arguments, [H HI] present 
model Hamiltonians produce only a small subset of the 
cuprate HTS properties a)-i). For example, the t — J 
Hamiltonian implementation of the RVB theory in the 
mean fleld(MF) approximation does not directly produce 
Tc, nor does it give the correct SC gap over the doping 
range of the SC state, nor does it give the observed ratio 
of the maximum pseudogap to the SC gap at optimal 
doping. [34] 

The Bardeen-Cooper-Schrieffer(BCS) phonon medi- 
ated electron-pairing theory of superconductivitv [451 146] , 
in conjunction with the method of Eliashberg. [47| ap- 
plicable to the strong interaction regime, [15] provides 
the framework for understanding the microscopic inter- 
actions responsible for LTS. However, the BCS phonon 
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mediated pairing interaction does not produce cuprate 
HTS, although the BCS theory appears to apply to 
MgB2 with Tc = 40°K.|49] Multiple pairing interactions 
are considered necessary to explain cuprate HTS,[H [13 
Uni mi |3| but the exact nature of the doping dependent 
interactions in the cuprate unit cell presently remains 
beyond quantitative measurement. Hence, formulation 
of a HTS model Hamiltonian relies on reasonable con- 
jecture about the underlying mechanisms responsible for 
HTS, with subsequent validation requiring comparison 
of model predictions with many diverse experiments. 

Our Hamiltonian H{x) is based on a phonon- 
mediated interaction V{x), detailed in Section HI, and 
an exchange interaction U(x), with the doping depen- 
dence X determined by the SPM. The U(x) < is pro- 
portional an effective spin-singlet exchange, which is re- 
lated to the X dependent particle occupation probabil- 
ity of the 02p states. The diagonal matrix elements of 
U produce significant static fluctuation(deviation from 
the MF) even for relatively weak U. We show that 
this static fluctuation is the key ingredient for HTS. For 
weak U, the model predicts the hsted a)-h) cuprate dop- 
ing dependent characteristics with a second order phase 
transition(SOPT). The SPM coupled with H produces 
a unique relation between the conjoined model and the 
observed phase transition boundary Tc{x), the SC gap 
A(r, x) and the NS pseudogap Apg{x). For stronger 
U, large static fluctuation results in an SC state that 
is essentially controlled by the ratio U/V. If U/V ^ 1 
the model exhibits a first order phase transition(FOPT), 
and quantum criticality, a phenomenon of broad interest 
beyond cuprate HTS systems. (SUl [SI] 

2 STATE PROBABILTIY 
MODEL 

Intrinsic cuprates are antiferromagnetic insulators with 
single or multiple CuO xy-planes, alternating with car- 
rier reservoir planes which may have a significant effect 
on the value of Tc{xop) in various cuprates. [24l l52] The 
Cu3d (Cu^"*") ion is in an octahedral environment sur- 
rounded by six 02p (0^~) ions, with the apical oxygens 
on the z-axis. Jahn- Teller distortion along the z-axis 
lowers the energy of the system by increasing the Cu 



apical O distance. The resulting intrinsic cuprate is a 
quasi-2D antiferromagnetic insulator, exhibiting strong 
electron correlation. In the a;y-plane nearest neighbor 
Cu3d and 02p orbitals point directly toward each other, 
producing strong covalent bond coupling. This inter- 
mediary ligand coupling indirectly connects neighboring 
Cu3d ions, giving rise to a large antiferromagnetic su- 
perexchange interaction Jm 0.13eV. 

A salient property of cuprate HTS is its dependence 
on the doping concentration x. For many cuprates, hole 
doping produces a dome shaped SC phase boundary in 
the doping range [xi — 0.05, a;2 = 0.27] with Tc{xop ~ 
0.16) exceeding 100°K; whereas for electron doping the 
range is comparatively narrow, [xi = 0.14, 0:2 < 0.2], 
with much lower Tc{xop)- This section introduces a 
general state probability model based on particle occu- 
pation of the cuprate electronic states. The SPM de- 
fines the probability of the SC state, and various nor- 
mal states, in terms of the probabilities of the accessible 
unit cell(UC) states, independent of the model Hamilto- 
nian. The SPM mandates the existence plaquettes with 
net charge and spin to preserve local charge-spin neu- 
trality. Relating particle occupation probabilities to the 
doping concentration x, the SC gap |A(x)|, the pseudo- 
gap Apg(x), and a new anti-pseudogap IS'pgix), which 
determine domain boundaries of the phase diagram, are 
a natural consequence of the SPM. Linking the SPM to 
our Hamiltonian determines the doping dependence of 
the thermodynamic quantities. Electron doping, with 
prediction of a possible new class of HTS materials, is 
considered at the end of this section. 

Although there is considerable hybridization of the 
atomic states, it is advantageous to denote each UC 
state by the occupation of the constituent single particle 
states. This simplification is consistent with a funda- 
mental assumption of quantum mechanics that compos- 
ite systems retain the properties of the individual con- 
stituent particles to a considerable extent. [55] 

The undoped cuprate state tpAF is charge neutral, 
Cu^"*"— with a hole on the Cu3d-orbitals and two 
electrons filling each of the 02 p^; and Pj, orbitals. For 
hole doping, the AF state, with few exceptions, [S3] is 
rapidly quenched by doping induced hopping. [51 In a 
doping range < a; < a;i there is an onset of charge 
transfer excitations involving random hopping of parti- 
cles among Cu3d orbitals and 02p orbitals. Initially 
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holes are localized, but as doping increases hopping pro- 
duces numerous possible states with net charge and spin. 

The electronic state of a doped cuprate is modeled 
here by UC states containing a total of 5 fermion par- 
ticles with spin a = ±1/2. These single-particle states 
form a large set of 2^ particle states x2^ spin states 
= l}^ = 1024 states denoted as (/?. It is advantageous 
to divide the UC states into two sets. The sets ^pha and 
^Pe(J represent states with a hole or an electron, respec- 
tively, in a given Cu3d orbital. The UC states, written 
as elements of a matrix, are 



(1) 



(fihaiij) = \ha)[<^><^>'^]^j, 



\ea 



, ha"), \ ha\ 



\ha' , ha"), \ ea' , ea") 



The single particle Cu3d states are denoted by \ea) and 
\ha). Elements of the column vector $, and the row 
vector <i>-^ are states representing the possible particle 
configurations for the p^-orbital, and the pj.-orbital, re- 
spectively. The matrix 4><I>^ gives the complete set of 
particle configurations, for each spin set. The tpha states 
are illustrated in table [ij 

Contrasting the well defined undoped AF state, the 
doped state is a probabilistic mixture of the states in 
Eq. ([T]) that result from doping induced hopping. The 
constituent single particle states are assumed to be mu- 
tually exclusive, and thus the distinct UC states are mu- 
tually exclusive, analogous to a 16 sided die with each 
face imprinted with one of the particle configurations, 
including spin. Let Pea-ix,T) represent the energy av- 
eraged probability that a given single particle state is 
occupied by an electron with spin a, at temperature T, 
and doping concentration x. The corresponding hole oc- 
cupation probability is Pha = 1 — Pea-- The probability 
that a particular UC state fhaiij) exists is a joint prob- 
ability involving products of the probabilities Pem and 
Pha- With one high T exception in the normal state, 
the scaled properties developed below have negligible T 
dependence. 

The states fha-i^j) and ^Pea{ij) carry net charge-spin. 
To maintain charge-spin neutrality in the CuO plane, it 
is energetically favorable for plaquettes containing UC 
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Table 1: Illustrated are the distinct single-particle states 
comprising the unit cell states iphaiij)- Holes are de- 
noted by o, and electrons by •, with spin not shown. 
The corner position represents one of the Cu3d orbitals. 
The undoped, charge neutral AF state is Other 
states have net charge, e.g. (p{12) has charge -t-e, (p{22) 
has +2e, and v(24) has -|-3e. The UC states ipeaiij) are 
obtained by interchanging every • and o, with the same 
spin. The net charge of ipeaiij) and ifha{ii) differ. 

states to form with opposite charge and spin for each 
value of X. (Plaquette is used in a generic sense to in- 
clude various geometrical shapes commensurate with lat- 
tice symmetry.) The sets of plaquette states necessary 
to preserve charge-spin neutrality, denoted as 



Lph-\, Ve-\, Vhl, fei, 



(2) 



must satisfy the state probability equalities 

PiiPht) = PM, PM = (3) 

where P{(p) is the probability that ip exists. Having es- 
tablished the UC states and the probability requirement 
Eq. ^ for pairing the plaquettes, one can determine SC 
and normal state probabilities and average properties. 

SC State: In a doping range xi < x < X2 b, coher- 
ent doping dependent SC state ipse emerges for T < Tc- 
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An effective exchange interaction responsible for the SC 
state is formulated below. Based on repulsive energy 
considerations, it is assumed that the SC state excludes 
the AF state t/5(ll) and the extreme overdoped state 
(^(44), with completely filled or completely empty p- 
orbitals, respectively. In the overdoped range x ^ X2 
states form in the charge transfer gap, resulting in a 
Fermi-liquid state (ppL with negligible remaining indi- 
vidual particle character, and the SC state is destroyed. 

We define the probability that the SC state exists as 
the sum of the probabilities of the UC states (p/icr that 
remain after extracting the high energy states V'(ll) and 
(/3(44) from Table [l] and the corresponding ip^a states. 
The result is 



PeaFp, Psci^ha) ^ Ph.Fp (4) 



the resultant HTS properties are clearly not the same, as 
discussed below.) The doped holes fill oxygen orbitals in 
the reservoir planes, which in turn distort the lattice and 
induce particle hopping in the CuO conducting planes. 
Since the probability of the SC state, Eq. ([5|, reduced to 
the that of effective p orbital states with an electron and 
a hole, it is reasonable to assume that the probability 
Ph that a position labeled h has a hole is Ph oc Ilh{x). 
Without loss of essential information, we set the con- 
stant of proportionality to unity, giving the probability 
Pe = Ile{x) = 1 — Ilh{x) that a position labeled e has an 
electron. Hence, Eq. ^ written in terms of the dopant 
probabilities is 



Pscix) = meix)Unix)[i - -n,{x)nh{x)], (6) 



(Pea' + Pea"){Pha' + Pha"] 



2P,„,P,,„Pha'Ph 



Adding the expressions in Eq. Q gives the SC state 
probability Pgc = Fpi which depends only on the occu- 
pation probabilities for the p-orbitals. For the remainder 
of the analysis of the SC state we neglect spin depen- 
dence of the p-orbital probabilities Pea' and Pha' , giving 



PSC ^ ^PePh[l ~ ^PePh]. 



(5) 



The first term in Psc arises from a simplified three par- 
ticle state model by considering the row and columns 
in table [T] independently, excluding doubly occupied p- 
orbitals. The second term results from a joint probabilty 
involving both and pj,-orbitals. Thus Psd as de- 
fined above, automatically excludes all doubly occupied 
p-orbitals. Writing Psc in the form Psc ~ ^ ^ Pe ^ Pit 
confirms that ip{ll) and f{4:A) are excluded from the SC 
state. 

The utility of Psc is implemented by relating the 
particle occupation probabilities to the doping concen- 
tration X. Assuming a uniform probability density over 
a doping range w — X2 — Xi, the probability that a 
doped hole is created in doping range [xi,xi + x] is 
Ilh{x) = {x — xi)/w. (Setting dli/dx constant, and ne- 
glecting temperature dependence, is the same assump- 
tion used in Anderson Ref. [34], but its application and 



^h{x) 



Xl 



Ueix) 



X2 



X2-X1, 



with Iie{x) + lVh(x) = 1. 

Equations Q-Q are independent of the detailed in- 
teraction responsible for the SC state. The proposed 
pairing interaction is a net effective inter-plaquette, spin 
singlet exchange J(x) per unit cell that emerges as a 
consequence of doping induced hopping. The J{x) is an 
average over the values for each UC state Lpij participat- 
ing in the SC state. Determination of the exchange Jij 
between states iph^(ij) and (feiiij) requires microscopic 
analysis to determine the overlap integrals and the re- 
sultant eigenstates of a very complicated system. Here, 
the relative value J{x)/ J{xop) is obtained by probabil- 
ity arguments, with J{xop) a parameter found by fitting 
experimental data in Section V. 

The exchange energy, per unit cell, averaged over the 
accessible UC state configurations ipij is formulated in 
Appendix A. The result from Eq. (63 1 is 



J(X) « J2li-e{x)Iih{x)[l - 2{\ - ^)llellh] (7) 

ZJ2 



Jl — >/ll + >/22 + 2( J12 + J34), 



J2 



2(Jl3 + J23), J3 = 2(Ji4 + J24), 



The Jij are the exchange constants corresponding to the 
interaction between iph^{i,j) and ipei{i,j). Thus J{x) 
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is an average spin-singlet exchange between the paired 
plaquettes that ensure local charge-spin neutrality. Al- 
though microscopic evaluation of the Jij involves multi- 
ple particle-particle and particle-antiparticle interactions 
between UC states, the doping dependence of J{x) is 
via the n(a;)'s which refer to p orbital occupation. The 
participation of the Cu3d particle serves as a mediator, 
appearing only in the exchange constants. Note that 
for J1/J2 = 3/2, the effective exchange is identical to 
PsCj to within a constant of proportionality. Since the 
maximum of {l/2)Ile{x)Ilh{x) = 1/8 the small Px-Py 
interaction correction will be neglected below. In this 
approximation the Jij in the factor J2 have no effect on 
scaled HTS properties. 

Retaining the linear two particle terms in Eqs. ([6]) 
and ([7]), the effective exchange J{x) and the SC state 
probabilities are given by 



gives the doping dependent relations 



(8) 



The simplification to the dominant two particle term im- 
plies that an effective electron-hole pair within the p^^ 
and Pj, orbitals is an essential ingredient characterizing 
the SC state. Assuming a lower energy alternating spin 
state, the implied quasi-particles are spin-excitons with 
spin-singlet exchange. This picture is somewhat distinct 
from the Zhang- Rice [35] and Geballe[24] pictures which 
explicitly contain a Cu3d particle. All three scenarios 
involve dynamical processes, which are treated here as 
static, time-averaged, phenomena. 

In view of Eq. Q, we propose an hypothesis: 

Doping dependence of scaled energy parame- 
ters that characterize SC, and normal, states 
is manifested only via the doping dependent 
probabilities that the relevant UC states are 
accessible. 

The validity of this hypothesis, already evident in Eq. 
([s]), is further substantiated by the following analysis 
that produces the doping dependence of many observed 
cuprate properties. 

In accordance with the hypothesis, the scaled SC gap 
A(a;) cx: Ilh{x)Ile{x) . This result is also confirmed by 
combining Eq. ( [46| , which is a direct consequence our 
model Hamiltonian developed below, with Eq. (Isl). This 



Tcjx) 
Tc{xop) 



A{xl 
Ao 



AIlh{x)Ile{x), Ao = AKp). (9) 



These universal relations, independent of the average 
density of states Nq and the cutoff temperature Tm de- 
fined in Section IV, give the doping dependent phase 
boundary Tc{x) and the low temperature SC gap A{x) 
observed in cuprates, as shown in Section V. Equation 
([9]), concomitant with the Hamiltonian model relation 
Tciiy), is supported by the (T, x) dependence of the Hall- 
coefRcient used to track the hole (electron) charge char- 
acteristic throughout the SC phase. [25] The two- particle 
nature of the pairing interaction parameter v{x) and the 
resulting A (a;), is consistent with measurement of the 
SC gap, requiring two-particle probes. [7] In the RVB- 
Hubbard model, Ref. [34], the SC gap is assumed to 
be proportional to gt{x) = [2a;/(l -I- x)]^, where gt is 
a kinetic energy renormalization factor. This function 
only approximates A (a:) in the very underdoped region 
X « Xop, whereas Eq. ^ agrees with the observed gap 
over the entire doping range of the SC state. 

Normal State: The doping dependent normal state 
is characterized by unusual properties, e.g. pseudogap, 
vortices, stripes, etc, reviewed in Ref. IT:. The apparent 
complexity is daunting, as are the myriad of complex and 
exotic theories. Nichtsdestoweniger, it is our contention 
that there is a rather simple explanation for many prop- 
erties based on the UC states in Eq. ([T]) and illustrated 
in Table dl 

Since doping induced hopping is a random process, 
it is asserted that the normal state is characterized by 
the complete set of UC states in Eq. ( pj). Assuming that 
all of the state configurations in Eq. M are accessible, 
the probability P{(ph) that some one of the states in ip^ 
exists is denoted by 



Pi^ha) = = PhaP{^)P{'^>^)- 



(10) 



The probability that <I> exists is defined as the sum of 
the probabilities for each state, i.e. it is the probability 
that some one of the states exist. Thus 

F($) = (Pea' + Pha'){Pea" + Pha") = 1- (H) 

Since P($)-'" is also unity, one obtains 

P{^ha)^Pha{x,T), P(^e.) =Pe.(x,T), (12) 
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independent of the p-orbital occupation for any given 
spin set. One can derive the same result for P{(ph,cr) by 
tediously summing the probabilities for each of the UC 
states Lpfic, in Table [T] and for the corresponding ip^^ 
states. 

The seemingly trivial expression in Eq. (12), a con- 



sequence of the completeness of the UC p-orbital states 
used, has profound implications. It gives a non-zero 
probability for static charge and spin fluctuation on any 
given Cu3d orbital. Since Ph(x) ^ Pe{x) for x ^ Xop, it 
is evident that the states (ph and cannot exist alone. 
As concluded above, plaquettes form in the CuO plane 
with opposite charge and spin for each value of x. The 
presence of such plaquettes is consistent with the for- 
mation of charge and/or spin density waves with con- 
comitant gaps. P [55] The paired plaquettes dictate that 
the system is composed of the four states in Eq. ([2]). 
Since Eq. ( [l2| shows that the doping dependence of 
these states is characterized by the Cu3d orbital occu- 
pation, we set P{x) — n(x), giving the plaquette state 
probabilities 



(13) 



P{ip 



hi) 



ne(x), p(^e;) = n,,(a;) 



Consider a checkerboard pattern of UC plaquette states 
tp^, and ip\^. Since ne(a;) -|- Hh{x) = 1, the plaquette 
state ip^{x < Xop) is dominated by spin-up electrons, 
and the plaquette state y^(x < Xop) is dominated by 
spin-down holes. At optimal doping x — Xop there is no 
net charge or spin for the plaquette pair. Overdoping 
gives the reverse of the underdoped picture. 

Since the probability of the normal states reduced 
to simple, one particle probabilities for the Cu3d or- 
bitals, we are able to formulate the exchange interaction 
between plaquettes in the normal states more precisely 
than that in the SC state. Let Jq = J{h t, e |) = - J(e f 
, h 4,) denote the spin singlet exchange between spins 
in the plaquette states characterized by the correspond- 
ing particle and spin in the Cu3d orbital. Neglecting 
particle-particle exchange, expected to be much weaker 
than particle-antiparticle exchange, and noting that ex- 
change requires joint probabilities, the doping dependent 
average interaction energy Js{x) for the balanced plaque- 



tte pair is 

J,(x) = M^lix) - nl{x)] = Jo[n„(a;) - n,(x)]. (14) 
Introducing the definitions 



l^h{x) = Jo + Js{x) = 2jQXlh{x) 



K{x) = Jo- J,(2;) = 2Jone(x), 



(15) 



the energy Ae(x) becomes identical to the experimen- 
tal pseudogap Apg(x) with the single requirement that 
Ae(a;) > A(a;). This gives the maximum SC gap Ag = 
0.5 Jo such that Ae(a;) is excluded from the SC state, 
yielding the form 



Apg(a;) = Ae(a;) =4Aone(x). 



(16) 



The pseudogap Apg{x) in Eq. (|T6| decreases linearly 



with doping from a maximum Apg(xi) = 4Ao to 
Apg(x2) = A (2:2) = 0. It is shown in Section V that both 



Eqs. ([9]) and ([Te]) are in excellent agreement with a broad 
class of cuprates. In comparison, the RVB-Hubbard the- 
ory in Rcf. [33], Figure 2, also gives Apg{x) > A{x), but 
Apg{xi) > 6A0 does not agree with experiment. 

A simple SC gap-pseudogap relation is obtained by 
eliminating the scaling factor Aq from Eqs. ^ and ( 16 1, 
or alternatively retaining Ao and writing a difference re- 
lation. This gives 

A{x) = Apg{x)Uhix) = Apg{x) - 4Aonl{x), (17) 

which is a universal relation corresponding to Eq. (j9]). It 
is evident that Apg and A have the same symmetry, and 
that their difference is proportional to the probability for 
filled p orbitals, which characterize the AF state. 

It is useful to introduce the concept of a pseudogap 
state, and an anti-pseudogap state, defined by their prob- 
abilities Ppg{x) = He (a;), and Pp„{x) = n;i(a;), respec- 



tively, corresponding to the pseudogap A^ 



Ae{x) 



and anti-pseudogap Apg{x) — A^(a;). The transition do- 
mains between the pseudogap(anti-pseudogap) state and 
the SC state is defined by the joint probabilities 

Ptrix) = Ppg{x)Psc{x), PUx) = P'pg{x)Psc{x) 

(18) 

where Psc{x) is given by Eq. ([6]). The corresponding 
energy transition boundaries are Atr{x) = AoPtr(a;) and 
A',,{x) = AoPU^). 
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Figure 1: Plotted is the doping concentration x depen- 
dence of the common, two-particle SC gap Ae/i(x) = 
A(a:;) and the distinct, single-particle pseudogaps Ae(x) 
and A;i(a;) for hole doped cuprates and the proposed 
electron doped anti-cuprates, respectively. The dashed 
curves are the transition energies ^tr{,x) and IS!^^(x). 



The SPM exhibits reflection symmetry about for 
energies in both the SC and normal states. The symme- 
try is illustrated in Figjl] where the SC two-particle gap 
A(a:) = Ae/i(x) and the distinct, single-particle, pseudo 
gaps Ae(a;) and A/i(a;) are plotted in units of Ag ver- 
sus the doping concentration x. The dashed curves are 
Atj.(a;) and Aj^(a;), representing the boundary of the 
transition from the normal state to the SC state. Not 
shown, but of interest, is the net normal state pseudogap 
energy i?pg(a::) = (l/4Ao)[Ag(a;)-|-A^(a;)]. As a; increases 
on the underdoped side Epg{x) and the pseudogap de- 
crease and the SC state becomes more robust. The SC 
state is maximized at optimal doping where Epg{x) is 
minimum, with equal pseudogaps and no net plaquette- 
plaquette interaction energy in the normal state. On the 
overdoped side Epg{x) and the anti-pseudogap increases 
and the SC state loses coherence. 

A pertinent question is why is the pseudogap Ae(x) 
observed, while Ah{x) has not yet been observed. In the 
overdoped range, it is plausible that a local excess hole 
imbalance would be obscured, or destroyed, by the onset 
of the fermi- fluid hole state, but Ah{x) should appear 



close to the SC gap on the underdoped side near Xop- In 
any case, Ah{x) has a significant role in both cuprates 
and the anti-cuprates proposed below. 

The doping dependence of the Knight-shift, which is 
proportional to the spin susceptibility Xs, is partially ex- 
plained by Ae{x,T). Here the width of the probability 
distribution w(T) — X2{T) — xi{T) is assumed to be a 
function of T, with Xop = {xi + X2)/2 constant. Since 
Ae{x,T) is proportional to the spin-singlet exchange be- 
tween plaquettes, Xsix,T) cc l/Ag(a;,T), giving 



= 1 



w{T) 



(19) 



In the underdoped range Xs{x,T) < Xs{xopiT), and in 
the overdoped range Xs{x, T) > Xsixop, T) for each value 
of T as observed in La2-i:Sra;Cu02.[SS] The T depen- 
dence of u'(T) is not obtained within the SPM. However, 
using a Curie- Weiss form w{T) = T — O approximates 
the Xs{x,T) curves. 

For hole doping of cuprates, a few more comparisons 
are noteworthy. The slave-boson picture developed from 
the t — J model, [2] divides the phase diagram into do- 
mains by a fermi pairing parameter and a MF bose con- 
densation parameter, which have the same doping de- 
pendence as Ae(a;) and A^(a;), respectively. The Nernst 
domain is very similar to our pseudogap-SC state transi- 
tion domain, bounded by Atr(x) and the SC gap A(x). 

In the Emery-Kivelson model [44] Xop is determined 
by the intersection of a pairing amplitude Ao(a;), which 
decreases with x, and the phase stiffness Eg(x) = ksTg, 
which increases with x. The value of Tc{xop) is gener- 
ally not the maximum value of Tr- fa:). [24j In our model, 
the intersection of the pseudogap Ae(a:) and the anti- 
pseudogap Aii{x) occurs at Xop, with Tc{xop) the maxi- 
mum value. 

Measurements of the superfluid density (T, x) [S] 
provide another connection to our model. The x de- 
pendence ns{T = 0,x) cx X — xi is the same as 
the that of anti-pseudo gap A'pg{x) oc Uii{x). Also, 
dns{T,x)/dT cx Hg(a;), which is the probability that a 
p orbital is full, fits the data as well as the singular func- 
tion dugiTjx) / dT cx x~'^ proposed.^ 

Electron Doping and Anti-Cuprates: Electron 
doping of cuprates produces an SC state in the compar- 
atively narrow range, [xi = 0.14,^2 < 0.2], with much 
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lower Tc{xop) than for hole doping. Referring to the 
discussion above, the SC state of cuprates stems from 
the states in Table [T] excluding ip{ll) and <p(44). Hole 
doping increases the probability that these states exist at 
the expense of the definite, undoped AF state In 
contrast, while electron doping may induce hopping, it 
also serves to maintain the predominance of the undoped 
state to higher values of x, thus suppressing the forma- 
tion of the states necessary to form the SC state. At a 
doping induced hopping level which finally produces an 
SC state, the onset of an electron Fermi-fluid state ifpL 
again quickly suppresses the SC state. 

Electron doping of cuprates is not a very effective 
means of producing HTS. However, envision a possible 
class of " anti-cuprates" with the same crystal structure 
as a cuprate, but with exchanged anion and cation roles: 
The Cu^"*" and the are replaced by ions A^~ and 
B^+, respectively. The SC gap is the same as that for 
cuprates, but the role of the pseudogap is now taken 
by the anti-pseudogap Aii{x). If anti-cuprate material 
exists, or can be synthesized, electron doping should ex- 
hibit HTS properties similar to hole doping of cuprates. 
The symmetry of the two HTS systems is evident in 

Figg 

3 Hamiltonian H 

The Hamiltonian is constructed from a combined inter- 
action V + U. As analyzed in Section II, a consequence 
of particle hopping is an induced effective spin-singlet 
exchange interaction U{x) oc J(x). However, it is em- 
phasized that the development in this section, and many 
of the resulting HTS properties, are completely inde- 
pendent of the X dependence of U . In addition to J7, 
the formation of small polarons contribute to a phonon- 
mediated interaction V . The resulting two-particle in- 
teraction Hint has the general form 



Hi, 



Pkq = 



(20) 



q kk' 



-k+q/2iCk+q/2tj 



where the interaction matrix T = [F^k'] = [Vkt' + ^^kk']- 
In the cuprate SC state it is consistent with ob- 



planes. [571 ISH] Accordingly, Hint is approximated by re- 
placing all terms in the operator products in Eq. (20) 
with BCS pairing terms, i.e. neglecting q dependence. 
The Hamiltonian for the electron system, including the 
kinetic energy is 



H = H. 



kin 



(21) 



Hkin = X! ["^k'^k + C^k^-k] , 
k 

H^nt{T) = P^rp ^rkk'(c_kCk)^C_k'Ck', 
kk' 

where p = [c_kCk] is a column vector with BCS pairing 
operators as the elements. The fixed spins are implic- 
itly indicated by ±k = k t, — k |. The single parti- 
cle energies, referenced to the chemical potential fi, are 
Efc = Ck — In the underdoped regime ek is the tight 
binding kinetic energy, and in the overdoped regime it 
is expected that Ck is characterized by an effective mass 
m* . The general form of T is retained in the analysis 
below until it becomes necessary to develop the specific 
matrix structure relevant to cuprates. Next, a method 
is introduced for dealing with non-negligible static fluc- 
tuation due to Fkk 7^ 0. 

Mean Field and Static Fluctuation: Our pro- 
cedure begins with an extraction of the mean fleld(MF) 
part Hmf of H leaving a deviation(static fluctuation) 
Hd- Then Hmf is diagonalized and H^ is approximated 
by (Hd), which is subsequently evaluated exactly in the 
eigenstates of H„^f. The first step is to introduce a vec- 
tor pair deviation operator d = [dk] = b — p, where 
d = [dk] and the components of b = [6k] are complex 
scalar fields 6k (x). Using d to reorganize HintiX) in Eq. 
(pll) gives the form 



H = Hk,n + Hb{T)+HdiT), 



(22) 



Hh 



Hd = 



E 

kk' 



Ekk'Kc-k'Ck' + (c-kCk)^&k' - K^k' 



kk' 



^kk'rffcC^k' 



If Hd is neglected, H reduces to the form of the conven- 



servation that BCS pairs c_k^Ckt form in the CuO tional BCS MF Hamiltonian with bir = (c_i^| ci^^). [581155 
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The MF approximation is valid for interactions with 
only off-diagonal matrix elements, such as the often used 
model with Fkk' = — Vo(l — <5kk')- ttiis case, static 
fluctuations are zero in the random phase approxima- 
tion and is negligible. However, as shown below, Hd 
is not negligible for a pairing interaction with diagonal 
elements Fkk 7^ 0. This key feature of F produces HTS. 



The Hamiltonian in Eq. ( 22 ) is a reformulation of H 
in Eq. (21), with no approximation. The Bogoliubov- 
Valatin canonical transformation [SSJ [5D] diagonalizes 
Hmf = Hkin + Hh{T) in terms of quasi-particle num- 



ber operators Tik = 7k7k> and gk = Aj^Ak- To complete 
the diagonalization of H, the deviation term is ap- 
proximated by its ensemble average (Hd) — (d'''rd). The 
details are outlined in Appendix A. Exact eval uati on of 
(dj^dk') in the eigenstates of H,nf, using Eqs. (851 and 
( 88 1 , gives the diagonal Hamiltonian 



H 



Hmf + {Hd) 



(23) 



k 

{Hd) = 5]rkk'[(4>('^k'> + J^kk'fT^], 



CTk = 



1 - 1^ tanh(/?£;k/2), p = l/iksT), 



£k + |AkP. 



The quasi-particle excitation energy E]^ depends on the 
energy gap Ak, which is linked to the average gap devi- 
ation (i5k)- They are defined by 



k' 



-X!^kk'(dk' 



(24) 



Using (dk) = ^k — (c-kCk), and Eq. (81) to evaluate 
-kCk) , the gap Ak and gap deviation (Jk) in Eq. ( 24 ) 



are related by the constraint 

ffk = l + ^5Irkk'(c-k'Ck'), (25) 



Ak 



(c-kCk> = ;^^tanh(^£;k/2). 



Setting git: = 0, Eq. ( 25 ) reduces to the conventional BCS 
constraint that determines Ak. However, it is shown be- 
low that Ak cannot be determined from Eq. ( 25 ) because 
((5k) and Ak are in phase and increase simultaneously 
for an interaction with diagonal matrix elements. Since 
< (c_kCk) < 1/2, it is evident from Eq. (25) that a 



large gap solution has a corresponding large deviation. 
Using H in Eq. (23) gives the model expressions 



for the thermodynamic functions defined in Appendix B. 
The thermodynamic potential and the average internal 
energy U — {H) are 



mf 



(26) 



mf 



^{- ln[2 cosh(/?i?k/2)] - ek - Ak^^, 



and 



U = -J2[Ek tanh(/3^k/2) - £k - Akfo^ + {Hd) . (27) 

k 

It should be noted that since the evaluation of {Hd) in 
Eq. ([23| is exact, the U{T) is exact for aU T. Ahhough 
n(r) is approximate, it is expected to be accurate for 
T <$:Tc since fi(0) = U{0). The third thermodynamic 
function of interest is the entropy S(T). Inserting U and 
ft from Eq. ([26]) into Eq. (|92| gives 



S = 



^^{(2/;9)ln[2cosh(/3£;k/2)]- 



£:ktanh(,3£'k/2)}. 



(28) 



The standard Fermion gas form for S applied in LTS.[55] 
is obtained using the identity x = ln[/(— a:)//(a:)]. The 
error inherent in il accounts for the absence of the direct 
effect of {Hd) in S. Such dependence is indirect, via Ak 
determined self-consistently from the fixed point of ft. 

Random Phase Approximation: To resolve the 
system for non-negligible static fluctuation, we apply a 
random phase approximation(RPA) to {Hd) in Eq. (23 1. 



Assume the phases of the complex fluctuation compo- 
nents (dk) are random. The only contribution to {Hd) 
in an ensemble average over the phases is from the phase 
independent k = k' terms. The fluctuation terms in 
{Hd) due to the off-diagonal elements of Fkk' are zero. 
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The average RPA value of the complex ((5k) is also zero, 
but I (5k) I ^ 0. Noting that (4) = (dk)*, the phase 
average RPA expressions are 



LB, 



1(4 



k 

Ell 



1(4 



(4')|= 



(29) 



(30) 



It should be noted that {fid) can make a significant con- 
tribution to H if Tkk 7^ 0, even when the |(dk)P terms 
are negligible. 

Interaction Matrix Structure: It is the goal here 
to formulate a minimal structure for F = V + U which 
reflects the essential physics of cuprates. The U = [C/kk'] 
is formulated in Appendix A from a generic spin-singlet 
exchange interaction, with the effective doping depen- 
dent exchange U[x) cx J(x) given in Section II. 

Since the doping range of the SC state of most 
cuprates is relatively small (0.05 ^ a; ^ 0.27), it is ex- 
pected that a phonon mediated interaction V has a sig- 
nificant contribution due to the formation of tight bind- 
ing small polarons. Small polaron formation in cuprates 
has three important consequences: [55] 1) The lattice de- 
formation in response to the charge variation tracks the 
electron hopping between Cu3d and 02p orbitals. 2) 
There is exponential reduction in the electronic band- 
width. 3) Electron hopping requires emission or absorp- 
tion of phonons. In the same order, relevance to our 
model is three fold: 1) It is assumed that the symmetry 
of V is the same Cu3d-02p bond symmetry as that of 
the exchange interaction. 2) There is a cut-off for the 
kinetic energy. 3) The diagonal elements Vkk = 0. 

An additional consideration is the repulsive electron- 
electron interaction Vc included vclV — ^Vf^. Since 

< and 14 > 0, the net phonon interaction V{x) — 
— [|V^(x)| — Vc(x)\. As the material changes with dop- 
ing from an insulator to a metal, the |T4(x)| is expected 
to decrease with increasing x as the strong small po- 
laron interaction changes to a weaker Fermi-fluid-phonon 
form. [14] The Vc{x) is also expected to decrease with in- 
creasing x as screening increases. Since x dependence of 
the difference |V^(2:)| — Vc(x) is reduced, V is approx- 
imated by a constant. Although V involves both CuO 
and reservoir planes, 3D effects [52] are incorporated here 
only in the strength of V . 



In accordance with the above discussion, the interac- 
tion matrix elements are given by Eq. ( 78 ) , which is 



kk' 



-ro(x)(l-,5kk')V^k^k' -2C/o(x)5, 



kk', 



(31) 



ro(a:) = V^^U^ix) C/o(x) =4J(x). 



Gap symmetry is corroborated by numerous experi- 
ments on cuprates indicating a mixed s- and d2,2_j^2-wave 
gap, [21 [m E] and by general gauge and time-reversal 
symmetry breaking arguments. [3] As shown after Eq. 
( [75[ ) the symmetry factor i/jk can be d-wave, or s-wave, 
or a complex linear combination. All three choices have 
the same effect on the thermodynamic functions, which 
depend on IV'kP- However, a d-wave gap Ak — A?/'k, 
with -0^ = co^(kx) — cos{ky) enforces the Hubbard dou- 



ble occupancy restriction {cl^cl^, 
from the expression for (c_kCk) in Eq- ( 25 ) 



— 0, which follows 



Using Eq. (31 1 in the definitions in Eq. (24) gives 



Ak « AV'k, A = roJ2'Ptbk 



(4) ~ S^Ji,, 6^roJ2K{dk) 



(32) 



The approximate Ak and (5k assume that the sum of 
terms in A and S is much larger than the single term 
arising from the diagonal elements of F. Using the ap- 



proximate forms of Ak and (5k, in the constraint Eq. ( 25 ) 



gives 



± 



lAI 



E 



tanh(/3£;k/2). (33) 



Since g is real, only magnitudes \S\ and |A| appear in the 
relative fiuctuation ratio. It is shown that .g > for an 
SC state to exist. 

Using Ak ~ AV'k from Eq. (32 1 to eliminate the 
fok from Hmf in Eq. (23), and using Eq. (31) in the 



phase average RPA Eqs. (29) and (30), the Hamiltonian 
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assumes the form 



(34) 



Ho 

{Ha) 
where 
El = 



= - y^[£'k(i - »-k - gk) " £k] 



To ' 



-2(l + a)C/oi-L_[/gSi, 



^ 



E,i = ^iv'knwp, E,2 = E(i-i'/'kp)i(dk)p 

k k 

The remainder of this article is based on the Hamil 
Ionian in Eq. (34), coupled to the constraint Eq. (33) 



The general symmetry dependent problem is outlined in 
Appendix D, where it is shown that the model not only 
produces small Uq HTS solutions with the |(5p term ne- 
glected, but also extreme HTS solutions exist for weak 
interactions in the large fluctuation limit. It is also con- 
cluded that the symmetry factor ■0k does not fundamen- 
tally change the HTS thermodynamic properties which 
depend only on IV'kP- 



4 AVERAGED 
ANALYSIS 



SYMMETRY 



In accordance with the above discussion, the sums are 
now transformed to integrals, with replaced by the 
average (|V'kl)<"' = which sets a = 0. The interac- 
tions V and U may have different energy scales. How- 
ever, underdoped cuprates are narrow band insulators, 
with further narrowing due to the formation of small 
polarons. Thus electronic kinetic energies arc limited to 
values much less than the polaron cut-off at T = 0. In 
the overdoped domain the band widths are larger but 
the effective pairing interaction is reduced, resulting in 
a small value of Tc- Since the use of multiple cut-offs 
would not essentially change the HTS results over the rel- 
atively narrow doping dependent range of the SC state. 



we invoke a kinetic energy cut-off. This simplifying as- 
sumption is implemented by replacing electron energies 
£k by their k-space angular averages (ek) , which are then 
bandwidth limited by (sk) < £m = ksTm- 

In this section we focus on the resolution of the sys- 
tem when the effect of the fluctuation term with |5p in 



Eq. (34) can be neglected. However, since the gap |A| 



is determined from the thermodynamic potential fixed 
point equation, dfl/d\A\ — 0, the validity of the small 
fluctuation solution can only be assessed by retaining 
the |(5p in fi, and then neglecting it in the fixed point 
equation. 

The integral forms of the thermodynamic potential 
il in Eq. (26), the internal energy U in Eq. (27), with 



Hj) from Eq. (34), and the entropy S in Eq. (28) are 



an(t,(t)) = 
aU{t,(j)) = 
a(H,) = 



l-AtIo{t,^)+'^+a{Hd), (35) 
7 



l-2l3{t,c^) + ^+a{Hd), (36) 
7 

-2i^(^V--/i(t,0), 
7 / X 



5(i,0) - 2fcBx[2/o(i,0)-(lA)^3(i, </))]. (37) 
The constraint Eq. (|33|) assumes the form 

(38) 



±K=g{t,cl,) = l~jI{t,< 



The integrals /„(i, 0) are defined in Appendix E. The 
scaled temperature and gap are 



t 



T 



and the material parameters are 

7 = + t^^NqVo, ly ^ NoUo ^ ANqJ 

X = NoEm, a = l/(xem). 



(39) 



(40) 



The A'o is the average density of electron states in the 
energy integration interval [—£„,£,„]. The i] = jjp — r/c 
is the effective electron-phonon interaction parameter rjp 
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reduced by the repulsive coulomb parameter 77c • As dis- 



cussed before Eq. (31 1, there is at least partial canceling 
of the X dependence of 77. It is possible that other pa- 
rameters, cut-off energy e„, density of states Nq, and 
chemical potential /x are also functions of x. However, 
lacking explicit information, the model is kept reason- 
ably simple by neglecting doping dependence of the set 
[r], Em, Nq, n]. Excellent agreement with experiment in 
Section V confirms the validity of implementing this ap- 
proximation. 

Gap Equation: The thermodynamic potential ( 35 ) 
contains two unknowns, |A(t)| and \S{t)\. The con- 
straints g = 1 — 7/, and dil/d(ly^\t = determine g 
and ^, self-consistently, as functions of the parameters 



[r?, J/, x] at each value of t. From Eq. (351 one obtains 
the fixed point equation 
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7 



2iy 



7 



X 



0, (41) 



where the integrals I2 and /4 are defined in Appendix 
E. The first term is the mean field term, the second and 
third terms are from the \d\'^ fluctuation term, and the 
last term is due to the constant diagonal elements of the 
effective exchange interaction U. In the limit v = 0, Eq. 
^ gives g = 1- ril{t, 0) = 0, which is the LTS BCS 



constraint for non-zero t. 



Neglecting the fluctuation terms, Eq. (41 1 reduces to 



7(^,0) = -/2(t,0)>O, 



1 

7 ' ■ X 
which is valid for parameters satisfying 



g_ 

7 



2V -l2{t,(t))+h{t,(j)) < 1. 

X 



(42) 



(43) 



(42) exist for 7 
0. 



= ?7 + 1/ > 0, since 
Thus the exchange 
Vc < 0, 



Solutions of Eq. 
the integrals / > and I2 > 
interaction v is so effective that even ii rj ~ rjp 
i.e. when the coulomb interaction overrides the phonon 
interaction, a SC solution exists for values of v as long 
as -q -\- V > Q. 

It is elucidating to consider approximate implicit so- 
lutions of Eq. (42) for the interaction limit 7/[l— ,g(i)] < 
0.3 at i = and t = tc- Using Eqs. ((TTil) and gives 



(/)(0) = 2e~i/''e3(°'/'^, tc = LlMe'^'^ e^'^'"^'^'^ . (44) 



The BCS MF results for A(0) and Tc are obtained for 
g{v = 0) = 0. As increases from zero, the g{t) > 
initially results in an exponential 5/7 increase in A(0) 
and Tc- For larger v the dependence of A(i) and Tc on 
v is linear over a relatively broad range of v values. 
The linearity is found analytically under certain con- 



ditions. It follows from Eq. (42), using Eqs. (114) and 
Eq. ( [TT7| that 



g(O,0) 
7 

7 



1 - 7r/4 V 



for 0(0) < 1 



(45) 



X 



for tfntc 



where 0.0738 < Ci{t) ^ 0.2. It is worth noting that e,„ 
and A^o cancel in both expressions in Eq. ( 45 1 ; hence the 



enhancement of the gap and tc, relative to the BCS MF 
values, is due only to the exchange interaction cx C/q > 0, 
independent of the cutoff energy and the average density 
of states. Using g{0) and g{tc) from Eq. (45) in Eq. 



(44 1, leads to approximate linear relations for the ratio 
of the t = gap and tc at v and i>o- They are 



V 



(46) 



which are independent of Aq and Tm- The importance 
of this positive linearity to HTS, in stark contrast to 
the negative exponential dependence in LTS, cannot be 
over emphasized. It was shown in Section II to be an 
essential characteristic for the doping dependence of the 
gap |A(a;)|, and Tcix) observed in cuprates. It is also 
found numerically that Eq. ( 46 ) holds over a relatively 



broad range of including values that give significant 
fluctuations requiring the solution of Eq. (41). Using 
Eq. ([8|, and setting = v{xop = xi -\- w /2) = . 
doping dependence of the effective exchange energy is 



the 



v{x) 



AIVh{x)Iie{x). 



Thus |A(a::)|, and Tc{x) are linked to x. 
Critical Field, and Specific Heat: 

sation energy 

A17(t,0) = VL{t,(t)) - f^(t,0). 



(47) 



The conden- 



(48) 
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defines a thermodynamic critical magnetic field He by 
(l/2)Moi/2 (t) ^ |Ar!(t»|. (49) 



At low temperatures, T < Tc, for > Eq. ([99|) yields 
the linear temperature dependence 



HcjT) 
Hc{0) 



B{<t>) 



1 - 



T 



(50) 



V ln4-l Tc 



It is shown in section V that for v > 10^"^ the slope B{4>) 
lies in the range [0.4, 0.6] , in agreement with cuprate HTS 
observation. The v = Q LTS limit in Eq. ( [99| gives a 
very different Hc{T)/Hc{'S) = 1 - 1.06(r/T^. 

The entropy determines the specific heat C — 
TdS/dT for fi constant. Using the integral relations 
pl3| in Eq. ([37l gives the form 



y2_ 



t d(j)' 
2lk 



cosh 



(51) 



The normal state specific heat C„(t) = C{t, 0) is given by 
the first term in Eq. (51 ) with Y = y. The discontinuity 
AC{tc) = C{tc,4' -> - Cn{tc) at tc, obtained from 
Eq. $ll is 



AC(tc) - -fcsxtanh 



2tr 



dt 



For v > Q there is a large increase in the discontinuity 
due to the large slope d(lP /dt\t^ < 0. A quantitative 
comparison between the HTS and LTS C{t) is given in 
section V. General relationships between C and other 
forms of specific heat are developed in Appendix C. 



NUMERICAL SOLUTIONS 
AND COMPARISON WITH 
EXPERIMENT 



The numerical solution of Eq. ( 42 ) is facilitated by start- 
ing at t = 0, using the exact integrals in Eq. (114) and 



then increasing t. The minimum free energy SC gap A(t) 
is shown in Fig. [2] Contours are the condensation en- 
ergy An(i) in Eq. ([48| divided by x^m- Contours to the 




T/T 



Figure 2: SOFT solutions for (j){t) and g{t) = \5{t)/A{t)\ 
of Eq. (l57l) plotted for 77 = 0.25, v = 0.035, and x = 0-01. 



The contours are the scaled condensation energy. 

left (right) of the zero contour are negative(positive). Fig. 
[2] is plotted for a cutoff parameter x = 0.01, phonon pa- 
rameter r\ = 0.25, and a comparatively small exchange 
parameter v — 0.035 = 0.14?7. The significant points 
are 0(0) = A(0)/fcBT„, = 0.392, <c = TclT„, = 0.164, 
and gap ratio /S.{())/kBTc — 2.39. For comparison, us- 
ing the same values of x and 77, Fig. [3] is the v = 0.0 
(52) LTS MF solution, with (?!)„/(0) = 0.0366, ^cm/ = 0.021, 
and /S.{Q) / {ksTc) = 1.76. The enhancements relative 
to the LTS MF values are A(0)/A„/(0) = 10.7 and 
Tc/Tcrnf = 7.81. Thus a 12K LTS becomes a 94K HTS. 
The remarkable effectiveness of the diagonal matrix el- 
ements of the exchange interaction is clearly evident by 
comparing these enhancements with the modest mean 
field only enhancement Amf{T = 0,77 + v)//S.,nf{T = 
0,7]) = Tcmfir] + v)/Tcmf{ri) ~ 1-63. Large enhance- 
ment for all t < tc is concomitant with relatively large 
values of g{t) = \S{t)/ A{t)\, which varies from g{0) = 
0.52 to g{tc) = 0.45. 

The SOFT tc is plotted in Fig. |4] as a function 
of v, and using Eq. (47) as a typical cuprate second 
order phase transition boundary tc{x) for hole dop- 
ing, im [21 [7] An abrupt transition to the insula- 



tor(metal) state occurs for x just outside the [xi,X2] 
range. The tc{x) plotted in the doping range [0.05, 0.27] 
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Figure 3: LTS mean field solutions plotted for r/ 
0.25, = 0.0, and x = 0.01. 




for 77 = 0.25, = 0.14?7, where tcii^) is a linear 
function of v, has the parabolic shape observed in 
rather broad collection of cuprates: Bi2Sr2CaCu208+5, 
YBasCuaOy^a, TlaBasCuOe+i, HgBa2Cu04+5 which 
have same Tc{xop) ^ 90 — 95K, but with different num- 
bers, n = 1,2,3, of Cu-layers per unit cell.[7j For this 
value of w = X2 — xi = 0.22, the curve is often referred 
to as the empirical " universal curve" . A more appro- 
priate designation is "universal parabolic". The general 



(531 



form given by Eq. (|9]) and the equivalent Eq. 
is a fundamental cuprate characteristic based on state 
occupation probability arguments. Other cuprates ex- 
hibit a parabolic phase boundary, but with a different 
values for [xi,Xop,X2], and maximum Tc- For example, 
other measurements on the two-layer Bi2Sr2CaCu208+5 
give [a;i « 0.06, Xop ~ 0.12], [SI] but the shape of 
the SC phase boundary remains parabolic. The phase 
boundaries of the single layer Bi2Sr2LaCu06+5 and 
Bi2PbSr2LaCu06 + S are also parabolic with [a;i w 
0.11, X2 ~ 0.235]. t62j 

The lower Tc{x) for these cuprates with Tc{xop) < 
40°K is obtained in the model by reducing the value 
of 77 and/or v, and changing the range of x. However, 
the scaled curves are independent of 77 and are given by 
below. Doping dependence of iVo(x)[62] and 
may play a role in determining the SC 
and they probably intro- 
Also, 



Eq. 



(53) 



0.05 0.1 0.15 
V and X 



ti{x)\ 

phase transition range [xi, X2] 
duce some asymmetry into the tc{x) boundary, 
the sign change of dii{x)/dx contributes to the observed 
asymmetry of the phase boundary for hole versus elec- 
tron doping. 

The T — doping dependent SC gap A(x) and 
the pseudogap ^pg{x) are plotted in Fig. [5] for -q = 
0.25, — 0.14077- The gap A(a;) has the same shape 
as Tc(x), since both are linear functions of u{x). The x 
dependence of both A(a;) and Tc{x) is given byEq. (|9|, 
which can be written in the form 



Figure 4: The scaled transition temperature tciy) for 
soft's is plotted as a function of the interaction pa- 
rameter V, and using Eq.(47), as a phase boundary 



tc{x), with X the hole doping concentration in the range 
[0.05,0.27]. The curve is for rj = 0.25, x = 0.01, with 
Vm = 0.035. The maximum is Tc{xop)/Tm = 0.164. 



Tc{x) A(x,0) 



Tc{xop) A(a;op) 



w 



Xl) 



(53) 



where w — 



X2 



Xl. For the case xi = 0.05, w 



0.22, Eq. ([53]) is identical to the empirical. [Ml [551 E] 
doping dependent function observed in Bi2212, YB123, 
T12201, Hgl201, listed above, with Tcixop) ^ 90- 95K, 
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Figure 5: The SC gap A(T — 0,x) and the pseudogap 
Apg{x) are plotted in units of kBTm as a function of the 
doping concentration x for the same parameter values as 
in Fig. [2] The maximum gap is A{xop)/kBTm = 0.392. 
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Figure 6: Scaled thermodynamic critical fields Hc{T) 
are shown as a function of T/Tm for 77 — 0.25, x = 0.01 
and three values of v. The vertical scale is in units of 



0) « 21 ± Imev, and Ae^p{0)/kBTc ~ 2.5 ± 
0.15, as summarized in Ref. [?]■ 

The model A(0,a;) curve fits the experimen- 
tal Aexp{0,x) curve for a cutoff energy = 
Aexpixop)/(t)ixop) « 51mev, at Xop = 0.16 « 1/6. Using 
26mev ~ 300K gives T„ « 589K. The r,„ found 
Fig. lis T„ = Tciexp,Xop)/tcixop) « 92.5/0.164 = 
564K. The relative difference in T,„ is 4.3%, which 
is within the experimental error in Aexp(a;op,0) and 
the spread of Tc{exp). The effective density of states 
-^0 = x/^m ~ l/(5ev), which is that of a metal, i.e. 
Nq ~ l/ep- This is consistent with BCS supercon- 
ductor behavior at optimal doping. |2 |7] Using the Nq 
gives the exchange constant at optimal doping J{xop) = 
0.25Uo{xop) « 0.0438ev, and Vq « 1.25ev. Comparison 
with the undoped Cu3d-Cu3d exchange Jdd = 0.13ev 
gives J{xop)/Jdd ~ 1/3. Relating this result for the 
cuprates in Ref. [7] with the commonly used t-J model 
ratio Jdd/t ~ 1/3,(2] gives the geometric mean 



The pseudogap given by Eq. ( 16 1 is 



Jdd — \J tJ{Xop), 



for 



^op 



1/6. 



For these cuprates the competition between kinetic en- 
ergy per unit area and exchange energy is characterized 
by 2xopt « Jdd- 



Apg{x) = 4A(a;op) 



X2 



(55) 



Using the above A{xop) gives Apg = 42 ± 4mev, and 
_ Apg(xi) — 4A(a;oj,) = 84 ± 4 mev. This is within exper- 
imental error of the measured pseudogap Apg = 76 ± 4, 
extrapolated to xi — 0.05. [7] With this starting value, 
Apg{x) fits the data in the SC range 0.05 < a; < 0.27. 

The thermodynamic critical field Hc{T) is shown in 
Fig. [6] versus T/T^ and in Fig. [t] versus T/Tc- The 
curves m Fig. [§ with V = 0.07, are based on the full 
Eq. ( [ST] ). For T ^ Tc/3 there is significant linear slope 
for all ^{x) > 0, which has almost no dependence on x. 
The field He is related to the magnetic field penetration 
depth A(T) by a simple argument. In cuprates the exis- 
tence of vortices with effective penetration area ttA^ are 
pierced by constant quantized flux for T <C Tc- Thus 
the flux ttX^{T)Hc(T) is constant, and Eq. (50) gives 
,s the form 
(54) A2(0) , „,.^T 



= 1-5(0); 



(56) 



A2(r) ' '^'Tc' 

For the range 10^"^ « ^ « 0.035, we obtain a slope B « 
0.4 — 0.6. This is in quantitative agreement with the 
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Figure 7: Thermodynamic critical fields Hc{T)/Hc{0) 
are shown as a function of T /Tq for the same parameter 
values in Fig. [6j The v = curve is the LTS critical 
field of the BCS model, which lies slightly below that of 
the two-fluid model Hc{T)/Hc{^) = 1 - [T/Tcf. 



observed [Ml EH |68l ESI IZQl [TT] linear T dependence of 
A-2(r) with slopes B w 0.5 ± 0.1. 

In the SOFT domain the specific heat C{t) = 
tdS/dt increases with increasing h'{x). As a conse- 
quence the specific heat discontinuity at tc increases 
with hole doping, reaching a maximum at optimal 
doping, as observed. [5] The normalized HTS specific 
heat Cs(i/ = 0.035, </))/(2A:bx); and the normal state 
Cn{v, (j) = 0)/(2fcsx) are plotted in Fig. [s] The disconti- 
nuity at Tc is (Cs — C„)/C„ = 3.38. The corresponding 
discontinuity for the LTS specific heat C{i' = 0,4>) is 
1.43. Thus the model gives the observed, large compar- 
ative, difference between LTS and HTS specific heats. 

The HTS internal energy W(r, I/) shown in Fig. [9] for 
f = 0.035 exhibits anomalous behavior in contrast to the 
LTS internal energy U (T, 0) , which exhibits a monotonic 
increase with T for both the SC and normal states (See 
Fig. 3.3 in Tinkham|SH]). The unusual dip in U{T, v) in 
the range h « 0.07 < t < ^2 « 0.14 < tc = 0.164 in- 
dicates a temperature range [Ti,r2] of stronger effective 
exchange pairing than that for T < Ti. Entering the 
dip from the left, energy is transferred to non-electronic 
parts of the system. As T — >• Tc the system is absorbing 
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Figure 8: HTS specific heat C, calculated from Eq. ( [37|, 
is shown as a function of t for parameter values in Fig. [21 
The discontinuity at Tc is 2.36 times the corresponding 
LTS value. 
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Figure 9: The internal energy U from Eq. ( |36| ) is shown 
as a function of t for parameter values in Fig. [2] 
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energy as SC pairs are breaking. In the normal state 
T > Tc energy is again emitted, indicating some form of 
order stemming from the exchange parameter u. Satura- 
tion, due to the kinetic energy cut-off, occurs for T ^ Tq 
with U{T, p)->-l- v/x- 

The oxygen isotope effect on Tc is characterized by 
the coefficient aiso = —dhitc/dXuM, with M the oxy- 
gen isotope mass. Quantitative assessment of the effect 
may be comphcated and involve the effective carrier mass 
m* on the overdoped side.fOl [T5] However, making the 
simple assumption that dM oc —dx to preserve charge 
neutrality during hole doping, it follows from tc{x) plots 
that 

^iso is positive for x < Xop, negative for x > Xop^ 
and independent of x near Xop- Thus, without rigorous 
derivation, the model qualitatively gives the observed 
doping dependent evolution of the isotope effect on Tc- 

6 LARGE FLUCTUATIONS 

Comparisons in Section V to cuprate HTS are based on 
Eq. (42) in which fluctuations were neglected. Here we 



The gap expression is valid for (2/3)x < v < rj. It is 
evident that (j) is singular at v/rj = 1. The unique point is 
attained for v/x = Uo/sm > 2/3, and v/rj = Uq/Vq — 1, 
independent of the density of states Nq. Since % <C 1 
the point is reached even for extremely weak interactions 
- a QCP characteristic. Near the QCP, in the FOPT 
domain, both 5 and A are very large, and 5^1. There 
is no order parameter that approaches zero near the QCP 
at the FOPT, as is the case for a QCP at a SOPT.[32] 

In the QCP limit the scaled condensation energy 
aAr2(0, (t>) = ar2(0, (j)) saturates to the expression 



aAr2„ 



1 - 2 



1 



1 

2^ 



(59) 



consider solutions of Eq. (41 ) to determine the deviation 
from solutions of Eq. (42) for the case considered in 
Section V with the maximum v{xop) = 0.035, and to 



The corresponding critical field Hc{0) also saturates, 
and it follows from Eqs. (37) and ( 115[ ) that the en- 
tropy S{t) — )■ 0. Although it has been conjectured that 
the observed spike in Af2 near optimal doping 130] may 
signify proximity to a QCP, the comparisons in Section 
V indicate that cuprates are SOPT materials far from 
the large fluctuation regime that leads to a FOPT near 
the QCP. 

In the large gap limit the the large gap Eq. (58) is 
related algebraically to the three band t — J model. For 



show the effect of large fluctuations. Solving Eq. (41) ^ (2/3)x, Eq. (58) assumes the form 
for g leads to the key integral equation for (t>{t). It is 



Qit, 



1 



[l± 



(57) 



A3(0) 



Vn-Un 



t = vi^o, NoJ = •; , ; ;' . (eo) 



2 

2 hjt,^) 1 
xK{t,cj)) 2v 



Setting A = Ep — Ed, which is the cuprate charge trans- 



which is the symmetry independent equivalent of Eq. 
(1051 with a = 0. Eq. (57) determines g and self- 



consistently, as functions of [ry, x] each value of t. 
Large Gap Limit and t — J Model Relation: The 

general development in Appendix D shows the existence 
of a large gap solution in proximity to an interaction 
dependent singularity for any gap symmetry. Using the 
same Q <C 1 expansion of Eq. (57), and the i — > limits 
for the integrals given by Eq. (116), the (j>+{0) and (/(0) 
assume the asymptotic forms 



1] — u 



9 



1 



(58) 



fer gap, Eq. (60 1 is identical to the analogous relation in 
the three band t—J model for particle hopping between 
the Cu3d and the 02p orbitals.[2] The relation between 
the models is remarkable, but difficult to interpret, par- 
ticularly since the phonon interaction Vq is not in the 
t—J model. 

Large Fluctuation Numerical Solutions 
The solution (t)-{t) of Eq. (57 1 coincides with the 
solution 4>{t) of Eq. (|42| for small p = v/'q = Uq/Vq. 
For fixed [??, x]i increasing p increases (t'-{t) until the 
SOPT solution (j}-(tc) = is lost and a FOPT solution 
(/)+ emerges. The FOPT occurs with (p+itsw) > at 
the SN transition with maximum temperature tsN and 
gap (/)+(isAr) > 0. In the FOPT domain there is also a 
tc < tsN with (l)+{tc) = 0- For tc < t < tsN there 
are two solutions for each value of t. Hence, the FOPT 
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Figure 10: FOPT solutions (j)+{t) and \S{t)\/\A{t)\ of Eq. 
(Isfl are plotted for v = 0.96477, rj = 0.25 and x = 0.01. 



exhibits t dependent hysteresis with no applied magnetic 
field. In the limit p — >■ 1, both \il{0, 0)| and tsN saturate, 
but (f>{t < Isn) oo. Since this unique point exists for 
all t < tsN, and depends only on the ratio p = Uq/Vq, it 
satisfies the essential conditions of a QCP. _30j 

Fig. [To] shows a FOPT solution </)+(<) of Eq. (l57). 



Evident is the FOPT hysteresis, with the SN transition 
at 4>+{tsN) a. higher t than the NS transition at 
(j)-f-{tc) = 0. The maximum T saturates near T ^ 1 
as the QCP is approached. This gap-Tp decoupling con- 
trasts the behavior in the SOPT domain, where both the 
gap and Tc are increasing functions of v. For very large 
Tc the ratio Tc/Tp « x lO'^ - lO'^, consistent with 
low carrier densities in HTS cuprates, and 10^ — 10"* times 
the LTS Tc /Tp values. [3] The condensation energy con- 
tours near t — tend toward the QCP saturation value 
aAQo = —24.5 calculated from Eq. (59). 

SOPT phase boundaries are plotted in Fig. 11 The 
upper tcix), plotted for Um = 0.8477, lies just below the 
FOPT domain. The curve is slightly pinched near the 
maximum Tc{xop) /T^^ — 1.06, where Tciy) is nonlin- 
ear. At Xop the enhancement factor is 48 times the BCS 
MF value. For comparison, the lower curve is plotted for 
the same parameters used in the typical cuprate phase 
boundary in Fig. 4] The value of tc{xop) = 0.177 is 
about 7% higher than that in Fig. [4] where fluctuation 



Figure 11: Plotted are SOPT phase boundaries tc{x) 
versus the hole doping concentration x for rj — 0.25, 
X = 0.01. The upper boundary is for Um — 0.8477 which 
is very close to the FOPT region, and the lower boundary 
is for Vm = 0.1477. 

was neglected. A significant characteristic is the linear- 
ity of Tc{v) over a wide range of v values, which pro- 
duces the parabolic phase boundary Tc{x) observed in 
cuprates. 

The SC gap A(T = 0) is plotted in Fig. [I2]as A(i/) 
and as A(a;). Linearity of A(;/) over a relatively broad 
range of values, gives a parabolic A(a;). The upper 
curve A(a;) is for parameter values very close to the 
FOPT domain. As the QCP is approached, the increased 
nonlinearity of A(7^) gives the triangular shape, which in- 
creasingly fills the phase domain below the pseudo- and 
ant-pseudo gaps plotted in Fig. [T] The lower curve is 
for same parameter values used in Fig. [5j based on Eq. 
( [42| where fiuctuation is neglected. For the lower curve, 
A(xop) = 0.451 and A{xop)/kBTc = 2.55. The quanti- 
tative effect of the fluctuations is obtained by recalculat- 
ing the value of from the experimental data, as done 
in connection with Fig. [5] This gives = 537 from the 
experimental gap, and = 523 from the experimen- 
tal Tc- The relative difference in Tm is 2.6%, compared 
with 4.3%. The error in the gap ratio is 2% compared 
with 4%. Thus for ly = 0.035 neglect of the fiuctuation 
term is justified. 
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Figure 12: The SC gap A(r = 0, x) as a function of the 
hole doping concentration x, and as a function of the 
exchange parameter ly for same parameter values in Fig. 
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7 CONCLUSIONS 

The combined state probability-Hamiltonian model pre- 
dicts a rather large set of doping dependent HTS prop- 
erties listed in the introduction, and fundamental prop- 
erties of the normal state. The relative success is at- 
tributed to the following properties of the model: 

• The state probability model, defines the probabil- 
ity of the SC and normal states based particle oc- 
cupation of a large set of unit cell states. To pre- 
serve local charge neutrality in the CuO plane, the 
SPM mandates the formation of distinct charge- 
spin plaquettes. The SPM determines the doping 
dependence of the SC pairing interaction and the 
SC gap, and in the normal state it determines the 
pseudogap and an anti-pseudogap. 

• The procedure introduced to treat the deviation 
from the mean field Hamiltonian produces an ex- 
act expression for the internal energy U (T) = (H) . 
The thermodynamic potential ^(T) is exact at 
T = 0. 

• The phase average RPA value of the magnitude \6\ 
of static fluctuation from the mean field is not neg- 



ligible for any interaction with non-zero diagonal 
matrix elements. Values of \S\ and the SC gap |A| 
are determined self-consistently at the minimum 
free energy. 

• The spin-singlet exchange interaction U, with k 
independent diagonal matrix elements, is a very 
effective pairing glue, resulting in significant static 
fluctuation from the mean field state. 

• The SC gap and the critical temperature Tc are 
linear functions of the exchange interaction param- 
eter V (X U over a broad range of values of v, in- 
dependent of the effective density of states and the 
kinetic energy cutoff parameter. 

• Large static fluctuation produces quantum criti- 
cality with concomitant extreme HTS properties, 
although relevance to cuprates is tentative. 

The theory captures the key mechanism responsible 
for a large SC gap, high Tc, a large A(T = 0)/kBTc 
ratio, and a low temperature cx T. It also gives 
the doping dependent mechanism, coupled with H, that 
is responsible for the HTS phase transition boundary 
Tc{x), and gap A(x,0) for a broad range of cuprates 
exhibiting a SOPT with relatively high Tc- The doping 
dependent probability model ascribes a physical basis to 
a "universal parabolic" function for Tc{x)/Tc(xop) = 
A{0,x)/A{0,x op) satisfied by cuprates. It includes the 
empirical function[Bll [S51 [7] as a special case. The 
phonon interaction parameter rj oc V gives only a mean 
field contribution the HTS state. Universality stems en- 
tirely from the exchange parameter The form of 
I'ix) is further supported by the concomitant universal 
SC gap-pseudogap relation, which agrees with all ob- 
served parabolic phase boundaries. Fitting the broad 
spectrum of cuprate data in Ref. [7] further confirms the 
relevance and internal consistency of the theory. These 
non-trivial experimental signatures substantiate the es- 
sential role of U (x) as a pairing mechanism. Symmetry 
of the state probability model suggests the possibility of 
HTS in electron doped " anti-cuprates" , with the roles of 
the cation and anion reversed. 

A FOPT with T dependent hysteresis and large static 
fluctuation occurs in proximity to a QCP. The QCP is in- 
dependent of the gap symmetry function, the density of 
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states, and the energy cut-off parameter. The unique de- 
pendence on the interaction ratio Uq/Vq with the emer- 
gence of quantum criticality, may indicate the possibihty 
of extreme HTS with weak interactions. 

Although the theory predicts many cuprate HTS 
properties, a number of aspects are not included in the 
Hamiltonian. Doping dependence of parameters other 
than v{x), and explicit contributions from the charge 
reservoir layers are neglected. A fundamental evaluation 
of the state dependent contributions to the effective ex- 
change J{x), and the phonon mediated V{x) is lacking. 
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A INTERACTION FORMULA- 
TION 

This appendix contains two parts. The first relates to the 
exchange interaction, averaged over the UC states using 
the probability model in Section H. The second relates 
to the determination of the interaction matrix elements 
in the Hamiltonian in Section III. 

Average Exchange. Referring to Section II, 
we consider a checkerboard pattern of unit cell 
plaquettes(±) with opposite net charge-spin. The UC 
states in a given cell are denoted by ^pu^, ^Pe^i a-nd those 
in the neighbor cell are denoted by (fh\, '^ei^ which satisfy 
the probability constraints in Eq. ^ to maintain local 
charge-spin neutrality. Formulation of the exchange in- 
teraction between states in the paired cells requires an 
extremely complicated microscopic determination of the 
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overlap integrals for the 4 x 16 constituent UC states, 
which in turn determine the eigenstates of the inter- 
acting system. Here only the relative values of the ex- 
change for different groups of accessible states are formu- 
lated. In the context of the density functional approach 
with a Kohn-Sham exchange correlation potential. [55] 
the exchange values are assumed to be proportional to 
differences between characteristic energies Vhai'ij) and 
—Veaiij) corresponding to the UC states in Eq. ([2]), as 
illustrated in Table [T] In the SC state, excluding the 
UC states '/'(II) and v?(44), the average energies for the 
paired plaquettes are 



E. 



= p^^p,Ph[.hP,Ph + hPl + JiPfX (61) 



Jk 



vr 



= -P^PePh [J[PePh + J'2Pl + J'iPe] ,(62) 



V, 



[Wll +W22 +2(t;i2 +1^34)]'"' 



e, h 



2[wi3 -I- V2-i 



qa 



•^3 



2[wi4 + W24] 



qcr 



The designation ecr or ha in 14 is implicit in every Vij. 
The P'''s are Cu3d orbital occupation probabilities, and 
all P's are 02 p-orbital occupation probabilities. 

Adding Eqs. ^ and ([62|, defining J = Eh^: + Ee^; + 
Ehi + Ef.^, and noting that Jk = — J[. leads to the ex- 
pression 



J 



P.Ph 



= JiPePh + J2[Pe + Pi] - 



{J2-.h)[PtPl+PiPl] (63) 

Since Jk = Vjf^ — V^^ each energy difference 

Ji] = Vh^{ij) - Veiiij) = ~[ve^{ij) - Vhi{ij)]- (64) 

is interpreted as an effective spin-singlet exchange inter- 
action between the states iph^{ij) and (peiiij). Noting 
that J2 sa J3, and using (Pe + P,i)" 1, gives Eq. ([?]). 

Matrix Elements of U and V: The matrix ele- 
ments of [/ oc J as a spin-singlet exchange interaction 
is formulated. A singlet exchange interaction between 



spins at coordinates r and r' has the form 

H - - V2J , 



(65) 



Setting r' = r 4- R, the operator sj^^., creates a spin- 
singlet with a translationally invariant energy exchange 



constant Jr.r+R = J'R.- Eq- (65) and the transform of 



Hex to momentum space, with compaction to a single 
spin term is given without proof in Ref. |72| . We de- 
velop the transform details to point out an oversight in 
the determination of the diagonal elements, which are 
important for HTS. 

Defining the Fourier transform 



7=^CkiCTiexp[jki -r], i = l,...,4 

ki 



and the shift ki = -k + q/2, k2 = k + q/2, k3 = 
q/2, k4 = k' + q/2, the transform of Eq. (65) is 

— ^2Jkk'"kq'=kq 



=kq 



(66) 
-k' + 

(67) 



kk' 



-k+q/2;Ck+q/2-f 



-k+q/2tCk+q/2i 



'^kk' 



C-k+q/2iCk+q/2t + Ck+q/2iC_k+q/2t 
^jRCXpH(k-k')-R]. 



R 



The k space representation of H^x is simplified by writing 



J. 



kk' 



J, 



kk' 



"^kk" where 



•^kk' 



•^kk' 



-k,k' 



k.k' 



k,-k'i 



- 1° 
•^k,-k'- 



(68) 



Noting the second form of ^kq in Eq. ( 67 1 , it follows that 
terms involving J^^., do not contribute to Hex] whereas 
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terms involving J^, combine to give a single spin set form 



N 



(69) 



kk'q 



Pkq — C_k+q/2^Ck+q/2t- 



For a square lattice, applicable to the CuO plane of 
cuprates, we assume the nearest neighbor plaquette ex- 
change Jr = —J is an average value for the x and y 
directions. Setting — Ry = 1, the transformed ex- 
change parameter is 



kk' 



-4J[C, 



kk' -t- 'S'kk'] 



(70) 



kk' 



5, 



kk' 



cos cos k'^ + cos ky cos k'y , 



= sin kx sin k'^ + sin ky sin k'y 



Since Ckk + 5'kk = 2, the even part of Jkk' is 

^kk' = -4J [(1 - '5kk')Ckk' + 24k'] ■ 
Defining d and s-wave symmetry functions 

V'k — cos kx — cos fcj, , V'k — cos /Ca; -|- COS fcy , 

it follows that 



(71) 



(72) 



(73) 



This substitution of symmetry factors is used in Ref. [72] , 
however we note that the diagonal elements are not cor- 
rect. The contribution of 5kk to J^^^, was overlooked. 
This term, which produces the k independent diagonal 
element in Jkk', is essential to HTS. 

It is convenient to introduce a complex symmetry 
factor ■i/'k defined as 



V'k 



1 



V2 



[^k + #k]> 



(74) 



which gives the product 

V'k V'k' = Ckk' + *[cos ky cos k'x — cos kx cos k'y]. (75) 

It will be shown that all thermodynamic properties are 
functions only of [V'k P = Ckk- Furthermore, as stated in 



Ref. |72| . the s-wave and d-wave contributions to Hint 
are the same. 

Using the complex V'k, we replace the Jkk' with the 
complex exchange interaction 



C/kk' 



-C/o[(l-<5kk')V'kV'k'+25, 



kk'J 



(76) 



where Ua — 4J. As formulated in Section II, J is re- 
placed with an effective exchange J{x) that depends on 
the doping dependent particle occupation of the Cu3d- 
02p orbitals. 

In the phonon mediated interaction V, the main con- 
tribution is from tight binding, small polarons formed in 
response to the Cu3d-02p hopping of electrons (holes). 



Drawing from the discussion above Eq. (31): The sym- 



metry of Vkk' is the same as that of the exchange J, and 
the diagonal matrix elements Vkk = 0- The structure of 
the off-diagonal matrix elements, other than the symme- 
try factors, is a complicated function of k and k', which 
is modeled here by a constant —Vq. Accordingly, the ef- 
fective contribution of V to the electronic Hamiltonian 
is approximated by 



Vkk' = -1^0(1 -4k' )V^k V'k', 

jiving a, V + U interaction 



r 



kk' 



-(l/o + C/o)(l-4k')V'kV'k'-2C/o<5, 



kk' 



(77) 



(78) 



B DIAGONALIZATION OF H 

The Bogoliubov-Valatin canonical transformation, [59, 
I60| to a new set of Fermion operators 7k and Ak is 



7k = MfcCk - Wkclk, •^k = MkC-k + WkCk 



(79) 



Ck 



Mk7k + WkAj,., c\ ^ -v^-fi^ + UkXl, 



where the coefficients satisfy |ukP + |fkP = 1- The op- 
erators 7k(7k) and A|^(Ak) create(destroy) quasi-particle 
excitations consisting of a correlated electron-hole pair. 



Applying (79) to iJfem in Eq. (21), noting that £k 
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e_k , and using the anti-commutation rules, gives 



Hkin — ^ Ek [^k 



2uk^^k(Ak7k)'^ + 2(MkWk)*Ak7k] 



(|wkP - \vi,\^)inu + qu) + 2\vi,\ 



(80) 



where nu = 7k 7k, and gk — ^k^k are quasi-particle num- 
ber operators. To transform Hi, we use 



C-kCk = '"k'f^k(l - n-k - <7k) + 
«)^Ak7k - (wk)^(Ak7k)^ 



(81) 



givmg 



= y^{~(Mk«k^k + Mk^kAk)(l - rtk " gk) + 



[{vlYA^ - «)^AyAk7k - 

[(zik)'Ak - (z;k)'Ay (Ak7k)^ + 6^Ak}, (82) 

where Ak — — X^k' ^kk' ^k' ■ The off-diagonal terms Ak7k 
and (Ak7k)^ are eliminated from Hq =^ H/^in + Hi, using 
the coefficient constraint 

2ek"k^^k - ?^kAk + «kAk = 0. (83) 
The solution of Eqs. (831 and |wkp + |vkP = 1 leads 



to the relations 

(■Uk/wk)Ak = £kT-Ek, 2£;kUk'^k = TAk, 



(84) 



2^k|Mk|' - -EkT^k, 2£;k|vk|' =i^k±ek, 



where Ei, = y^el + \A^. Using Eqs. ((SSj) and (|84j) in 
Hkin + Hii{r) leads to the diagonal mean field form 

Hraf = J2^±Ekil - rik - Qk) + ek + Ak&y . (85) 



Determination of (Hd): To complete the diagonal- 
ization of H the operator dj^dk' is approximated by its 

average (dj^dk')- I^efine X[(pk9k')] = (Pk^k') - (Pk)(9k')- 
Using the definition dk = ^k ^ c_kCk and applying Eq. 
(ISll), the 6k 's cancel and one obtains 



-'^[((C-kCk)'^C_k'Ck')] 



x[(4dk' 



ui^vlul.,vi^,X[({ni^ + gk)(?^k' + Qk'))] 



(ukUk')^((Ak7k)^Ak'7k') + 



«i'k')'(Ak7k(Ak'7k')^) 



(86) 



The average of all other terms in X[{dl^di^r)] involving 
unmatched creation and annihilation operators are zero 
in the eigenstates of i?m/ ■ Using the anti-commutation 
relation for Fermion operators to rearrange the last two 
terms in Eq. (86 1, and noting that (7^/7k) = ^kk'("-k) 
and (Aj^,Ak) = 4k' (^k) gives 



X[(4dk/ 



WkWkUk'"k'^[((?^k + 9k)(n-k' + Qk'))] 



<5kk'|wk| ("k^k) 



<^kk'|wk|''((l - "-k)(l - gk)), 



(87) 



where (5kk' = for k 7^ k' and (5kk — 1- Since nk and (/k 



are uncorrelated it follows from Eq. (89 1 that (rik^k') = 
('^k)('7k') for all k and k'. Similarly, (nkrT-k') = (rik) {n-k') 
for k ^ k' , but (rik^k) = (nu) , and (f/k'Zk) — (gk), since 
the eigenvalues are and 1. Using these relations, one 
obtains ^[((rik + (Jk)("-k' + ^k'))] = 4k'[("k)(l - (^^k) + 



('7k)(l- (qk)]- Finally, noting that (?ik) = (gk), Eq. ([87]) 
reduces to 



X[(4rfk' 



r'5kk'(iVk)^ 



(88) 



It is stated in the literature that one should choose the 



lower sign in (84), but all thermodynamic functions are 



invariant with respect to the choice of sign, which is sim- 
ply a choice of an electron or a hole picture. 



^(^k) = |Mk|'(r^k) + K'k|'(l-(«k)), 

The (A'k) = (^k) ^^'^ average non-interacting fermion 
gas particle number density for state k for both spin 
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orientations. Using Eq. (|84|) in Eq. (88) and noting 



that (4rfk') = (^Xc^k') +^[(4^k'>] gives {Hd) in Eq, 



(23). In obtaining Eq. (23) wc apply Eq. (89) to obtain 
(«k) = (9k) = f{E^) = W^^ + l]-\ and use 1-2 f{x) = 
tanh(x/2). 



C THERMODYNAMIC FUNC- 
TIONS 

General relations between several thermodynamic func- 
tions are derived from their basic definitions, and differ- 
ent definitions of specific heats are related. 

Grand canonical ensemble average of an operator Q: 



(Q) = Tr{pQ), p^^eM-PH), 



Tr[exp(-/37?)], 



(89) 



where /3 = l/iksT), H{fi) = H{0) - pN with chemical 
potential piT), Z is the grand partition function, and p 
is the density operator. 

Thermodynamic potential (generalized free energy): 



= -(l//3)ln(Z). 
Von Neumann entropy; 

S = —kBTr{pln p). 



The above definitions (89) - (91) give the entropy 

s = -{{H)-n) = -— + {—). 



(90) 



(91) 



(92) 



The partition function for H in Eq. ( 23 ) is 
Z = [][l + exp(-/3£;k)]2exp[/3(£;k-ek)] 

k 

exp(-/3Ak6j;)exp(-/3(i/d)). 



(93) 



Using Z one obtains the expressions for Q in Eq. ( 26 ) 
and the internal energy U = (H) in Eq. (27). 



Specific Heat: Differentiating S in Eq. (92) with 



respect to T and eliminating S, yields the useful deriva- 
tive relation 

T^ = ^-(^) = -T^+tA(^) (94) 



The various forms of specific heat in Eq. ( 94 ) are 



C 



Cu = 



T 



d{H) 
dT 



= C 



'dT' 



(95) 



(96) 



(97) 



Details for C are given in Eq. (51 ). It is evident that for 
any model with (dH/dT) ^ 0, the specific heats differ. 
In this case the internal energy U cannot be determined 
from an integration of C with respect to T, as done for 
LTS.[5H] 

Condensation Energy at Low Temperature: 

The condensation energy is Ail{t, cj)) — il{t, (p)—fl{t, 0) < 
0, where n{t, 0) is given by Eq. (35) In the hmit 2t ^ 0^, 
and 2t ^ 1 for (f> = 0, the integrals in Appendix E give 
n(t, 4>) ~ r2(0, (p) to within exponentially small t depen- 
dence, and 

af7(i,0) - — (ln4- l)i- — (98) 
X 3 



Noting that AQ{0, 

An{t,(t)) 

A(0) = 



= ri(0, (j)), leads to the expression 
T / T 



1 - 2B{(f>) 



V In 4 



Tc 
1 Tc 



A{cj>) 



Tc 



(99) 



xWQ,4>)\T„. 



3|ar2(O,0)| 



T 



The factor in Eq. (99) shows that the linear T de- 



pendence is due to the diagonal matrix elements of U . 
For V — Q the LTS dependence (T/Tc Y' is recovered with 
A w (2/3)(7rTc/A)2 w 2.12. But, even for a smaU value 
of u the linear term in T/Tc dominates, since x ^ 1- 



Eq. (99) gives critical field Eq. (50) 



D GENERAL ANALYSIS 

The purpose of this appendix is to outline the behavior of 
the general, symmetry dependent minimum free energy 
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solutions. The corresponding to the model H in Eq. 

(100) 



(34) is 



n 



mf 



^mf + {Hd 

4 



{Ha) 



-eo^-C/oSi, Oo = 2(1 + a)(7o, 



(33) is 



where a ^ 1 is defined in Eq. ( 35 ) . The constraint Eq. 
±^=3=1- ToS, To = + C/o, (101) 



and the k-space sums are 

So = ^^ln[2cosh(/3£;k/2)], 

k 

^ - ^E^ta-h(W2), 



(102) 



^1 = jll'^k: 



ak = l-^tanh(/3£;k/2), 



where _Ek = a/ ^ w + l ^kp and Ak = Ai/jk- It follows 
from Eqs. ([l00|-(|l02|) that 







(103) 



P = 



a 9 In a 
1 + aain|A|2' 



The sums arising from the derivatives of Si and S, re 
spectively, are 



k 



|2 gfc 



k| -^^fkCk, 



2|Ak|% 
kl ^g2^''k, 



tanh(/3£;k/2) - ^ cosh-2(/3^^/2). 



The first term in Eq. ( 103 ) is from the mean field Hm/, 
and the remaining terms are from {Hd). Replacing the 
symmetry factor IV'kP in S^i, defined in Eq. (35), by its 
average value unity, a = and p = 0. Since a w 0, it is 
treated as a parameter and p ^ 1 is neglected. 
Solving ai7/a|A|2 = for .g/Eo gives 

1 ^ W \ 
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i± Vi-Q 



(105) 



W ^l- eoS4 



Solutions of Eq. (105) are the gap amplitudes |A±(T)|. 



Eq. ( 105 ) is complicated and it has several distinct so- 



lution domains depending on the relative values of the 
parameters. Extreme solutions of Eq. ( 105 ) are the 



small and large gap solutions that occur in the same 



limit Q ^ 1. To linear order in Q, Eq. (105) leads to 

To-eo 



Toeo 



= -(S-S4 



w' 



(106) 



E = 



(107) 



1 

Eq. (106) determines jA+j and Eq. (107) determines 
|A_|, corresponding to the sign in Eq. (105). Setting 



U{) = 0, Eq. (107) reduces to the BCS constraint for 
symmetry tpk, with a small, LTS gap |A_|. Since S — 
S4 > 0, Eq. ( 106 ) has no finite real solution in the limit 
Uo^O. 

Eq. ( 107), with W Ki 1, is the small Ua equation that 
follows from the MF part £7™/ plus the diagonal inter- 
action term proportional to Si in Eq. (103), neglecting 



terms generated by |(5p. As C/q is increased from zero Eq. 

has an effective T^ff = To/{l - f/oroSa > Eq that 
causes the gap amplitude |A_| to increase exponentially 



from the BCS value. [See Eq. (44).] The term small Uq 
is quantified by the condition 



eo[{/oS2 + S4] « 1, (108) 

required for negligible fiuctuation effect from |(5p. 

Retaining the fiuctuation terms leads to the emer- 
gence of a second solution |A+| determined from Eq. 
(104) (106), which becomes large when Fo — Oq ~ 0. In this 



limit the large gap solution of Eq. ( 106 ) is 

Toeo 



|A+(0)r = 



Vo-{l + 2a)Uo' 



(109) 
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(4/|A+P + V.2)3/2 



[UqCT]^ - £k] 



Existence of |A_|-| requires a minimum value of Uq such 
that > 0. Eq. (109) is singular when the interaction 
ratio p = Uq/Vq assumes the value pq = 1/(1 + 2a) ^ 1. 
Near the singularity the gap |A+(0)| — > oo, and the Em 
is essentially independent of |A+|. 

Although |A+| becomes increasingly large for p near 
Po, the g saturates to its maximum 5 = 1, and remains 
finite. It follows from Eq. (100 1, with some manipula- 



tion, that the asymptotic saturation value is 



(110) 



for any symmetry \ipk\. 

Several inferences are drawn from the analysis above: 
The symmetry factor ipi^ does not fundamentally change 
the thermodynamic properties. For small values of J/q, 
defined by pM] ), the |A_| is independent of \S\. The 
model Hamiltonian with |(5p neglected is applied in Sec- 
tion III, and its relevance to cuprates is clearly mani- 
fested by extensive comparison with experiment in Sec- 



tion V. When the inequality ( 108 1 is violated, retention 
of |5p produces a gap solution governed by p = Uq/Vq 
with a QCP at po ^ 1- Extreme HTS properties of the 
model are considered in Section VI. 



E INTEGRALS 



tion variable y = (ek) /em , are 
dy In 
t dio 



The integrals in Eqs. (351 - (38 1, and (57), with integra- 
tion vari 



hit, 



hit, 
hit. 



2 cosh 



= 2- 



1 



2 Jo 

1 dh 

(f) d(f) 

1 



dy 



1 — — tanh i — 



Y 



Y 

n) ' 

n 2 



= / dy 



Y 



Y2 



2t 



1 - — tanh ( — 
Y \ 2t 



(111) 



1 



— tanh — cosh — 



2t J 2t 



2t 



dyYi'dvii f — ) 

V2ty 



^ 







— tanh 




Y 





^ cosh ^ I ^ 
2t \ 2t 



Integration of I2 by parts gives the useful form 

2 

2/2(^,0) 



■ tanh 



2t 



dy 



The integrals satisfy the relations 







— tanh 


V2i). 


Y 





2t 



dh 
dt 



dt 



,dh 
' dt 



1 d£_ 

2 In" 



^ / dy 



Y' 



t 

2lH 



(112) 



(113) 
'Y' 



cosh — 
\2t. 
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For t <^tc the integrals are given by their i = hmit 
4i/o(O,0) 



7(0,0 



4(0, 



/i(O,0 



^2(0,0 



ri + 02/(0,0), yi = v/rT02, 



arcsinh(l/0) = ln[(l/0)(l + ^i)] 



= 1/n, 



(114) 



= l-Yi+0[l-(l/2)arctan(l/0)] 



-(l/ri)[l- 1/(2^1)] + 



(1/0)[1 - (l/2)arctan(l/0)]. 



For > 2i, 



2t/o(t,0) =/3(t, 



(115) 



Expansions of the integrals in Eq. (114 1 in powers of 
1/0 <C 1 are 



(116) 



/i(O,0) 
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2 20 
In the limit 0(t) ^ 

/i(i,0) = 1 -i[tanh(l/2i) +21ncosh(l/2i)], 

(117) 

h{t,0) = 0. 
For t < 0.1, with 7^ = 0.57726, giving D = 2.2677, 



/(t,0) 



InD 



tanh(x) ss ln(i:i/2t), 



/•C30 

/ (ixlnxcosh"^ X = In (4/7r) 
Jo 



(118) 



7e, 



Ci(t) 



1 1 f'/^'dx 



« C,{t) = --^-l^ ^tanh2(x). 
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